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ABS TRACT 

A method of estimating the depth of homogeneous mixing 
in the upper layer of a two-layer ocean is discussed. The 
method applies during the period of general heat loss from 
the ocean surface and is based upon convective processes 
occurring within the homogeneous layer. 

The model used is one developed by N. P. Bulgakov; 
density and salinity distribution with depth are plotted 
from the surface through the halocline. The density scale 
is oriented such that the maximum density for any given 
salinity is read directly. One can thus represent the 
maximum density obtainable (for a given salinity at 
freezing temperatures) by the salinity plot. By also plotting 
the initial density distribution with depth, and the 
average salinity attained in a homogeneous layer mixed to 
the plotted depth, one can determine the maximum depth 
of convective mixing from the intersection of these curves. 

The model was tested using data taken during the 
winter months over a five year period at Ocean Weather 
Station "P" located in the northeastern Pacific Ocean. 

The conclusions reached show that the model has two 
important features; first, it clearly demonstrates the 
important factors causing convective processes to occur; 
and second, it can be useful in separation of effects of 
the convective process from those of other mixing or non- 


mixing processes in the upper ocean layer. 
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‘ie Background. 

Knowledge of the vertical thermal structure of the 
oceans is a prime determinant of successful anti-submarine 
warfare. Maximum efficiency in the use of acoustic 
detection devices in the oceans can be achieved only 
through knowledge of its thermal structure. In particular, 
the ability to predict the depth of homogeneous mixing in 
the upper layer would greatly enhance the effectiveness 
of mobile sonar platforms. Toward this goal, investigation 
needs to be made concerning processes which help determine 
the ocean-thermal-structure distribution. In certain 
geographical areas during particular seasons the process 
of convective mixing is predominant in the upper layers of 
the ocean. 

This paper reports tests of a model for estimating 
limits of homogeneous mixing in the ocean due to convective 


mixing during the cooling season of the year. 





Za Introduction. 

The mixed-layer depth is defined as the depth to which 
the upper layer of the ocean has been mixed homogeneously. 
There are several causes for vertical mixing in the oceans. 
Basically they can be sub-divided into two categories 

1) mixing due to convective turbulence, and 

2) mixing due to frictional turbulence caused by wave 

action and currents. 
In addition there are two non-mixing types of fluctuations 
of the mixed-layer depth which are associated with 

1) internal waves, and 

2) the vertical convergence-divergence caused by 

changes in the overlying atmospheric pressure 
field, or by convergence-divergence in the surface 
current field. 
The possible influence of internal waves and convergence- 
divergence effects on the mixed-layer depth were not 
evaluated nor accounted for in this study. 

Convective mixing is due to vertical instability. In 
the cooling season there is a net loss of heat from the 
ocean surface caused by upward transfer of sensible and 
latent heat. The subsequent cooling of the water increases 
the density of the surface layer, thereby decreasing the 


vertical stability of the water column, or making it 





unstable. When instability is caused by cooling, convective 
mixing occurs down to a depth where the density is 
equivalent to the increased density of the surface layer. 
This process occurs until there is no longer a net radiative 
loss to the atmosphere or until the freezing temperature 

is reached. 

High rates of sensible heat loss to the atmosphere 
from the ocean surface usually are accompanied by high rates 
of evaporation. Evaporation acts to increase the surface 
density in two ways - by energy loss from the surface through 
the latent energy needed to cause the evaporation; and by 
increasing the surface salinity as the pure water is 
evaporated away, leaving behind in the ocean surface layer 
a greater concentration of dissolved solids. Therefore, 
convective mixing can develop from both sensible heat loss 
and a salinity increase at the surface of the ocean. 

The depth to which convective mixinz extends will 
depend to a degree on the initial density gradient present 
in the column of mixed-ocean water. When the initial 
density gradient is small, a given amount of cooling or 
evaporation will set off convection which will reach a 
relatively great depth. When the initial density gradient 
is large, convective mixing will extend downward a 


relatively small amount. For given cooling amounts, a 





certain critical density gradient can be determined which 
will not permit convective mixing to extend below a certain 
depth. 

The density after mixing of several ocean layers is 
greater than the weighted-average density of the layers 
mixed if a salinity gradient is present, due to the non- 
linear dependence of density on temperature and salinity 
(pressure effects can be ignored when compared to the 
temperature and salinity effects on density in the first 
100m of the ocean below the surface)(1}. This phenomenon 
has been called “contraction on mixing" by N. N. Zubov [2]. 

N. P. Bulgakov has proposed a model to determine the 
limiting depth of convective mixing for given initial 
conditions[3]. In this model he has attempted to account 
for this non-linear dependence of density on temperature 
and salinity, especially in brackish water areas in sub- 
polar and polar regions. The model gives a simple numerical 
or graphical determination of the maximum depth of 
convective mixing. According to this model the limiting 
depth of mixing without ice formation occurs when the layer 
has been cooled to its freezing temperature; any further 
increase in convective mixing might be caused by a salinity 
increase due to ice formation at the ocean surface. In the 


greatest part of the open oceans, the upper layer is not 
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cooled to its freezing temperature, and this limiting depth 
of convective mixing usually is not realized. 

Bulgakov's equation (to be discussed in detail in the 
next section) is the basis for this study which applies the 
convective model to data from a particular ocean region 
during months when convection is anticipated as a dominant 
mixing process. The model gives estimated depths of the 
homogeneous layer which may be verified against measured 
depths. Resolution of consistent differences provides an 
estimate of the validity of the model as well as a suggestion 
of the relative importance of convective mixing in that 


ocean situation. 
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2% Bulgakov's Model. 


In order to understand Bulgakov's approach to 


determination of the limiting depth of convective mixing, 


several definitions are needed. 


1) 


~ 
ba eal 


3) 


4) 


Density values are in terms of sigma values where 
sigma (© ) = (density - 1) x 103. 


Reduced density (O, ) is the maximum density for 
any given value of salinity (at the freezing 
temperature if the salinity is greater than 
24.6095 %o or at the temperature of maximum density 
if the salinity is less than Fis bO5 7 /aal\e. 
Convective mixing depth (H,. or H) which is the 
depth of the homogeneous layer after mixing due 

to convective processes (less than the depth to 
the bottom of the halocline). 

A one-layer ocean is defined as one where the upper 
boundary of the halocline (thermocline) reaches 
the surface, with consequent increasing salinity 
(decreasing temperature) commencing at the surface 


and continuing to a finite depth. 


A two-layer ocean is defined as one where the upper 
boundary of the halocline (thermocline) lies below 


a homogeneous mixed upper layer. 





5) Zero level (subscript O) which is the top of the 
halocline in a two-layer ocean; the ocean's 
surface in a one-layer ocean. 

6) Actual density (C,) which is the density 
corresponding to the temperature and salinity at 
any given level. 

Bulgakov's approach relates, for a vertical column of 
ocean extending through the halocline, the profiles of 
initial density distribution and average salinity after 
mixing (plotted against depth) in a basic equation{3]. 
Figures 1 and 2 are diagrams of the geometric representations 
for a one-layer and a two-layer ocean respectively. The 
initial salinity distribution is represented by the line 
labeled S; due to the corresponding GO scale, this line also 
represents the reduced density distribution, Or , prior to 
mixing. Assuming that salinity varies uniformly with depth, 
the line labeled S represents the mixed-layer salinity 
with changing depth of the convective layer as well as the 
reduced density of the layer after mixing; these values are 
plotted against the existing depth of convective mixing. 

The initial density distribution is represented by the line 
labeled Oy, . The important plots on these diagrams are the 
averaged-salinity curves and the initial density 


distribution curves. The limiting depth to which convective 
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processes can occur is determined by the intersection of 
these two curves. It is at this point that the actual 
density (in the initial density distribution) is equivalent 
to the reduced density (or maximum density) of the 
homogeneous layer after convective mixing. 

The geometrical construction of figures 1 and 2 provide 


a basis for the equation 


Se Os = B20 +S, (7) 


It should be noted that even though (1) is derived from the 
geometrical construction of the figures, each of the 
quantities has a physical interpretation. They are: 

1) On. = Oe is the difference prior to mixing 
between the density at the depth of convective 
mixing, Oy» and the density at the zero level; 

2) QO is the difference between the reduced density 
at the zero level, Or, , and the density at the 
zero level, O,, (surface or top of the thermocline) ; 
this term represents the temperature reduction 
necessary in order to mix the layer down to its 
maximum depth; 

3) A®s is the difference between the reduced density 
at the zero level, Or, » and the density at the 


depth of convective mixing, O\: 
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Derivation of the depth of conwmectivemmixime, Ho, fer a 
one-layer ocean based upon equation (1) follows. In figure 
l it is assumed that the initial salinity distribution 
varies uniformly with depth; therefore, the final salinity 


of the mixed-layer is 


S ast Sant Sie (2) 
gg 


Equation (1) can be transformed to 


Cyeemeoge 4p (3) 


The quantity Os is reduced to 


Cene GISS cs H ¢ 5 (4) 
where — is the salinity gradient in the halocline; 0.38 
is the change in reduced density per corresponding change 
in salinity. The quantity a8 represents the slope of the 
Cc 


dashed curve of figure 1. The constant, 0.38, is obtained 


by 
Orr cy” : 
Hea” Wo “a 
an — O,oOee 
Equation (1) may now be written as 
do = 
=- He = AG + 0.38 Ho (5) 


Itfhe constant 0.38 is determined by Bulgakov; this author, 
using Knudsen's Oceanographic Tables, determined the constant 
to be O.40 for the salinity range of 32.0 to 34.0 


aL 





from which the limiting depth of convective mixing, H can 


e? 
be evaluated if the temperature at the surface (to evaluate 
4&O), and the salinity and the density distribution are 
known. 
The derivation of an equation analogous to equation 
(5) for a two-layer ocean resembles that for a one-layer 
ocean, with a few important differences as figure 2 
illustrates. The derivation of a vertically-averaged 
salinity for the mixed-layer after convection is not as 
simple as in the one-layer case. Since Bulgakov does not 
give a detailed derivation for the two-layer ocean case, 
the following reasonable assumptions are made by this author 
in deriving the final equation stated by Bulgakov. The 
assumptions are 
1) the salinity distribution of the upper mixed layer 
is uniform with depth, 
2) the salinity varies uniformly with depth in the 
halocline, and 
3) the curve of the average salinity when plotted 
against maximum depth of convective mixing is 
essentially parabolic in shape within the halocline. 
Any weighted-average salinity for the convectively-mixed 
layer will give a much smaller value for the two-layer 


ocean case than for the one-layer case when convection 


tee 


‘ ‘ 
| | 

‘= <i eo —- . 
os» "1 _—_=— ee 


7 © oi © — « Men 
- a : : 
oe — —a » 3 et 
a oe ~~ 
= ee — + ' a, 
eo © oe eee ee ee — 


ie} to -« em es 


2 7 a « 


t—_ > = ————— h 





penetrates to a given depth into the halocline. The low- 
salinity water of the upper layer in the two-layer case 
acts to dilute the convectively-mixed layer. 

It can be shown that the two-layer ocean counterpart 


of equation (4) becomes 


Qo 
nN In 


ww 
Si eowos —— aw (6) 
3 


where H is the depth to the bottom of the halocline and 
dS 


aS is the initial salinity distribution with depth in the 


halocline. The equation comparable to (5) then becomes 


Clean ie dS H? 
=, H=AT + 0.19 _— (7) 


Equation (7) may be solved then for H, which represents the 
maximum change in depth of a homogeneous layer due to 
convective processes when all sensible heat above the 
freezing temperature has been removed from this layer. 
Equation (7) is the key equation in this study of convection 
penetration into the halocline; (7) will be referred to 
hereafter as either Bulgakov's model or Bulgakov's equation. 
Bulgakov's equation may be used as a tool for determining 
maximum changes in the depth to the top of the thermocline 
(assuming that the top of the thermocline corresponds to the 
top of the halocline). From equation (7) it is seen that 


H, the change in depth due to convection, depends on obtaining 


I 





1) the density and salinity gradients in the halocline, 

2) the actual and the reduced density at the top of 

the halocline, and 

3) the depth to the bottom of the halocline. 

The salinity gradients can be determined either from 
recent oceanographic stations or from climatology of the 
region. Actual density gradients may be approximated by 
using actual temperatures at the top and bottom of the 
halocline from bathythermograph (BT) observations and 
averaged salinities. Reduced density is evaluated initially 
as a function of salinity (which determines the freezing 
temperature} from tables; once determined, the reduced 
density scale may be plotted as an alternate scale for the 
salinity scale as in figures 1 and 2; thus the plot of 
salinity versus depth may be read also as the plot of 
reduced density versus depth. 

Thus equation (7) can be used readily to compute the 
maximum H (due to convection processes) for any area of the 
ocean where vertical salinity and temperature structure is 


known. 


1The temperature of maximum density was the freezing 
temperature due to the high salinity of the test area. Had 
the salinity been less than 24.695%., the temperature of 
maximum density would have been used. 
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lig Application. 

The northeastern Pacific Ocean was the area chosen to 
test the usefulness of Bulgakov's equation for estimeting 
the convective change in the depth of the mixed layer as 
described by the temperature structure (thermocline) in the 
region. The data were taken by the Canadian weather ships 
which occupy the ocean weather station "P", located at 
50° north latitude and 145° west longitude. 

The area has two distinct features which make it ideal 
for study. First, the station is located near the center 
of the Alaskan current gyral. Thus the advective change of 
vertical temperature structure due to permanent currents 
can be considered negligible. Secondly, the area is 
characterized by a pernanent halocline, located at a depth 
of about 100m at ship 'P'. The halocline separates the 
seasonally-affected upper layer from the nearly constant 
features of the lower layer',7. The upper mixed layer has 
an annual average salinity of 32.7%,,while the average 
salinity at the bottom of the halocline is 33.8%,[5]. 

The density structure of the area was considered to be 
a function of the teiuperature and salinity only. The 
adiabatic chenges were considered negligible as were the 
effects due to pressure change with depth, in accordance 
with Sverdrup's analysis of stability within the upper 100m 


of the ocean[l]. 


i) 





For a more detailed analysis of the physical structure 
of the northeastern Pacific, the reader is referred to 
Dodimead, Favorite and Hirano(5]}. 

The data supplied by the Pacific Oceanographic Group 
of Canada included both oceanographic station data and BT 
observations at ocean station "P"[7]. The data are 
considered to be extremely reliable. Eight oceanographic 
stations are taken (on the average) each month; and at least 
two BT observations are taken daily at O200z2 and 1700z hours. 

The data used were from November, 1957, through March, 
1962, for the months of November through March. Mean 
salinity gradients for each month were determined using the 
oceanographic stations (See values in Table I). The density 
gradients were determined for each observation using actual 
temperatures at the top and the bottom of the thermocline 
from BT observations and salinities of 32.7 and 33.8 % 
respectively. At the top of the thermocline the reduced 
density is 26.245 gm-cm-3. All densities are given as 


Sigma values, where 


Sigma = (densicy soles 
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TABLE 1 = Meayesalinzty Craduenes 


Month 2 ( Be) 
November OLA 
December rotten 
January = OES 
February 0268 
March OHS 7¢ 


The Control Data Corporation 1604 Computer was used 
along with the observed data mentioned above to solve 
equation (7) for the change in depth to the top of the 
thermocline, H, attributed to convective mixing. The 
density gradient, gs , reduced density, O. , and actual 
density, O, , were evaluated from the observed data while 
the salinity gradient, d= , was taken from the averaged 
data. After obtaining the quantity, H, this value was then 
added to the depth of the top of the thermocline to obtain 
the predicted depth of the mixed layer. 

As a second application of equation (7), the computer 
was used to solve for the density gradient using arbitrary 
values of H (lm increments), averaged values for AO and 
j° An analysis of various 
terms of equation (7) rewritten as 


2. , and arbitrary values for H 


dg = 2S dS 
ae. ~ -l = O.13 dl H; (Fay 
Ca) (2) 





; oe ommend be | ' sl scalie 
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shows that term (1) predominates in determining ae since 
the vertical salinity gradients are small in the open 
oceans as is the ratio 4, in tern (2). “Ateship “"Pgermn 
(2) is negligible in comparison with term (1) for the 


computation of cc , and equation (7a) may be written as 


dg . AD b 
= = (7b) 


Table II is based upon equation (7b) using lm increments of 

H and averaged values of AO (values are given in the table). 
This table, therefore, may be used to determine H by entering 
with a computed value of gs, The value of a may be 
computed from a current BT observation as outlined previously 
for application with equation (7). The interpretation which 
must be placed on the quantity ag is that it represents 

the critical value of the density gradient; at density 
gradient values greater than Sf , the convective processes 


will terminate at a depth equal to or less than the 


determined value of H. 
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os Interpretation of Data. 

Interpretation of data is highly subjective due mainly 
to the defining of a thermocline. This author interpreted 
the top of a thermocline as requiring at least a 0.5C 
decrease of temperature in 10m. The bottom was interpreted 
as the location where the temperature decreased less than 
eon On. 

A decision had to be made as to when only convective 
mixing is assumed to have occurred. Since the surface 
salinity changes at this station are very small, it can be 
assumed that evaporation effects on density were negligible. 
Convective mixing is assumed, thus, as a function only of 
heat loss in the surface layers during these cooling season 
months. From BT observations, it is extremely difficult to 
read the sea surface temperature to an accuracy greater 
than 0.1C. Therefore, it was assumed that when the sea 
surface temperature decreased by O0O.1C or more in any period 


of time, only convective mixing took place. 
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6. Determination and Evaluation of Results. 

Table III is a condensation of the statistical analysis 
of predicted mixed-layer depth in relation to observed mixed- 
layer depth. Approximately 1,100 BT observations were taken 
during the period of time covered in the analysis. Only 133 
of these observations were used, in the final analysis, as 
representative of purely convective changes. Bulgakov's 
equation for determining H in a two-layer ocean was used to 
derive the predicted mixed-layer depth. The predicted value 
was arrived at by adding to the mixed-layer depth at the 
time of prediction the computed depth change due to convection. 

Without exception, the arithmetic mean of the predicted 
mixed-layer depth, MLD,» was greater than the mean observed 
mixed-layer depth, MLD. Equation (7) predicts a maximum 
change associated with the removal from the layer of all 
sensible heat above freezing temperature; since this seldom 
occurs, a correction factor was developed to correct this 
over-prediction of convective mixing influence at the top 
of the thermocline. 

The over-prediction is not linear with temperature 
decrease in the surface layer which prevents the use of a 
constant correction factor. However, since the over-prediction 
is related to the retention of sensible heat, a ratio of 


change in sensible heat to initial total sensible heat was 


all 





TABLE III - Statistical Analysis of Data Results 


November December January February March 


ry 2323 .666 .678 pees fill 926 
ro 872 .668 698 1876 ~9L0 
MLD 68 100 108 digyzak 12h 
NID, 82 105 113 128 132 
Mj D. 68 LOu 106 120 slays 
SD eS T2e 1 Ome Pats Liles 
SD, 9.1 14.8 apes ae) qe 
Sil; 9.8 15 6.2 7.5 10.2 
IW 9.0 120 Ze) Iyings (9) SO 
r, =-  theslinear correlati onveom = pRedicted MEDeve 
observed MLD 
mo - the linear correlation of §unescorrected predicate: 
of MLD to observed MLD (meters) 
MLD ~ arithmetic mean of observed MLD (meters) 
MLD, - arithmetic mean of predicted MLD (meters) 
MLD, _ arithmetic mean of corrected MLD (meters) 
SD - standard deviation from MLD (meters) 
SD, - standard deviation from MLD, (meters) 
SD5 - standard deviation from MLD» (meters) 
LW - amplitude of internal waves at the mean observed 


MLD for waves of a period of 7 hours (meters) 


Le 





used to multiply the initially computed value of H. The 


corrected change in depth due to convective mixing is 


= ae 8 ~7,) (| 


Covrected Shes. com outed : 


H (8) 


i is the maximum annual temperature of the ocean surface 


max 
and T, 1s the temperature of the ocean surface at time of 
prediction. A climatological mean Tray of 13.40 was used 

in the correction factor [4]. Thus the correction factor 

is a funecion of 1, only. A correction facter Pomme ach i 
has been computed and is displayed in Table IV. 

The arithmetic mean of the corrected mixed-layer depth, 
MLD», was in close agreement with the observed mean mixed- 
layer depth, MLD. 

The linear correlation of predicted to observed mixed- 
layer depth both before and after correction are given in 
Table III. With the exception of November, the correlation 
coefficients increased as the cooling season progressed. 

In November, the great difference in the correlation 
coefficients before and after correction appears to involve 
the high temperatures of the surface waters during the 

early part of the month. An analysis of October was dropped 
because heating season effects were predominantly 
superimposed upon the convective or cooling season effects; 


the mixed-layer depth in October was too shallow to escape 


the frictional mixing of wind-induced turbulence. 
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TABLE IV - Initial Temperature vs. Corrective Factor 


LEout a weMe CommecGien Lnsateeasl Commestien 

Temperature, Factor Temperature Factor 

ikem tG.) Toe (e) 
Coen 694 8.2 388 
lin 686 Sr ~ 380 
Lee 679 Sik 25 
Led soe ee » 366 
te 664 Smo 358 
4.6 opi one 5 ge 
4.7 O49 8.8 348 
4.8 642 8.9 . 336 
LQ meswn ore) 226 
ore O27 SC e-wll we 
pret 619 0.2 ode 
OMe sole 9.3 BSWOle 
Dua 604 9.4 298 
vere SS i ie, Tage 
Dima> ~590 9.6 28h 
Fale) oer Shel WGN: 
Soe Be) 5 9.8 269 
Dro TOOy 9.9 208 
5.9 560 O50 iB 
6.0 De Orr 246 
aaa Die LOR 2 239 
One Boi ORS B25. 
O58 530 NONE weedy 
6.4 522 TOR Ay Als 
ont Sae5 ORG ~ 209 
6.6 507 Lol Ol 
Cn, 500 1.0 36 ee 
opus 492 10.9 er 
6.9 485 phi ve ny? 
fue 478 jai ie 
gill 1,70 Piles els 
ne 463 i li fall S 
73 Lp dl oe lie 
Wolo LL& Ties oie 
ac) LL0 rE ye mS 3 
(ede 433 de 27 
(aq E25 holes ~ 119 
Wats 418 11.9 re ona 
Leg 410 1220 ~ LOL 
SA 103 sb 097 
egal 396 LZ 32 089 
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The standard deviations were derived for the computed 
values of MLD, and MLD> to determine if they are of the 
same order of magnitude as the observed MLD; in Table III, 
it can be seen that the order of magnitudes are related with 
a maximum difference of 2.7m. 

The computed mixed-layer depth values, both uncorrected 
(MLD,) and corrected (MLD2), are shown by months and 
compared to observed values of mixed-layer depth (MLD) in 
figs. 3-7. The wide scatter of points in the figures tends 
to decrease as the cooling season progresses. Internal 
waves, it is believed, are the primary cause for the scatter 
of these data points. The amplitudes of the internal waves 
measured in the northeastern Atlantic by Glinskii and 
Boguslavskii range from 8m to 30m over a period range of 
7.0 hours to 68.0 hours; the amplitude of waves is a 
maximum at a depth range of 100-200m. As an example, the 
amplitudes at 80m for periods of 24.8, 7-0; 14.8, l12¥iy ama 
68.0 hours are approximately 8.0, 90, 9.5, 13 and 22m 
respectively(7]. Assuming that these results can be applied 
to a Similar region in the Pacific, the standard deviations 
for MLD, MLD,, and MLD5 in Table III match closely the 


amplitudes of the 7-hour and the 12.4-hour internal waves 
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Fig. 3'- Observed vs. predictéd/mixed-layer depth for November 
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F.lg.4 - Observed vs. predicted: mixed -layer depth for December 
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; Fig 5.- Observed vs. predicted mixed-layer depth for January 
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Fig G -Observed vs. predicted mixed=-layer depth for February 


29 





PREDICTED MIXED-LAYER DEPTH (METERS) 


100 10 120 130 140 150 1G0 170 180 .190 200 


90 





“COMPUTED VALUE ,MLD, 
CORRECTED VALUE, MLD, 


MARCH 
1958-1962 


1 
90 100 110 120 130 140 150 160 170 180 


OBSERVED MIXED-LAYER DEPTH (METERS) 


Fig. 7 - Observed vs. predicted mixed-layer depth for March 
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Ve Conclusions and Acknowledgements. 

It is apparent from the data results that Bulgakov's 
model can be useful in determining the maximum or limiting 
depth of convective mixing. 

The main value of his work is that it gives a suggestion, 
in terms of easily-measured parameters, of the mechanics of 
convective activity in the open ocean. Further study might 
be attempted to separate the effects of convective mixing 
and frictional mixing using Bulgakov's equation as a basic 
tool. Present work in the heating season tends to emphasize 
the wind-induced frictional mixing while ignoring the 
increased density (thus convective mixing) effects of mixing 
two or more layers of the oceans. Contraction on mixing of 
several layers should be accounted for in the heating season 
and Bulgakov's equation might be used since contraction is 
automatically accounted for. 

It can be concluded that Bulgakov's equation is useful 
in evaluating processes affecting the thermal structure of 
the oceans. 

I gratefully acknowledge the encouraging help and the 
critical comments of my faculty advisor, Associate 
rFrofessor G. H. Jung. My thanks to Dr. Tully of the Paci 
Oceanographic Group of Canada for the data made available 
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